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ABSTRACT 

A basic assumption in current halo occupation model is that the properties of a galaxy 
depend only on the mass of its dark matter halo. An important consequence of this 
is that the segregation of the galaxy population by large-scale environment is entirely 
due to the environmental dependence of the halo population. In this paper we use 
such a model to predict how the galaxy luminosity function depends on large-scale 
environment. The latter is represented by the density contrast (5) averaged over a 
spherical volume of radius R = 8 h~^Mpc. The model predicts that the Schechter 
function is a good approximation to the luminosity functions of galaxies brighter than 
~ 10^h~^ Lq (6j-band) in virtually all environments. The characteristic luminosity, 
L*, increases moderately with S. The faint-end slope, a, on the other hand, is quite 
independent of S. However, when splitting the galaxy population into early and late 
types, it is found that for late-types a is virtually constant, whereas for early-types a 
increases from —0.3 in underdense regions {S ~ —0.5) to —0.8 in highly overdense 
regions with S 10. The luminosity function at Li,^ < 10^h~^ Lq is significantly 
steeper than the extrapolation of the Schechter function that fits the brighter galaxies. 
This steepening is more significant for early-types and in low-density environments. 
The model also predicts that the luminosity density and mass density are closely 
correlated. The relation between the two is monotonic but highly non-linear. This 
suggests that one can use the luminosity density, averaged over a large volume, to 
rank the mass density. This, in turn, allows the environmental effects predicted here 
to be tested by observations. 

Key words: dark matter - large-scale structure of the universe - galaxies: haloes - 
methods: statistical 



1 INTRODUCTION 

In a hierarchical cosmogony like the cold dark matter 
(CDM) model, galaxies are assumed to form in dark mat- 
ter haloes (e.g. White & Rees 1978). A generic prediction of 
such cosmogony is that the properties of the galaxy popula- 
tion must be closely related to that of the halo population. If 
the cosmological density field is Gaussian and if the power 
of density perturbations extends to large scales, as is the 
case in the current CDM cosmogony, halo properties are ex- 
pected to be correlated to some degree with the large-scale 
structure, and so the properties of the galaxy population 
are also expected to change with large-scale environment. 
With high-resolution numerical simulations much has been 
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learned about the halo population. It turns out that a strong 
correlation with the large-scale environment exists only for 
halo mass; all other important halo properties are at most 
weakly correlated with the large-scale environment at any 
given time (e.g. Lemson & Kauffmann, 1999). This result, 
combined with our current model of galaxy formation, im- 
plies that the environmental dependence of the galaxy popu- 
lation is mainly due to the change of the halo mass function 
with large-scale environment. We can therefore hope to un- 
derstand the environmental dependence of the galaxy pop- 
ulation by understanding how galaxies occupy dark haloes 
of different masses. This is indeed the approach taken by 
the so-called halo occupation models, in which the number 
of galaxies per halo is assumed to depend on halo mass only 
(e.g. Jing, Mo & Borner 1998; Peacock & Smith 2000; Sel- 
jak 2000; Scoccimarro et al. 2001; Benson 2001; Bullock, 
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Wechsler & Sommerville 2002; Berlind & Weinberg 2002; 
Kang et al. 2002; Zheng et al. 2002; Yang, Mo & van den 
Bosch 2003a; van den Bosch, Yang & Mo 2003; Maghoc- 
chetti & Porciani 2003; Scranton 2002; Zehavi et al. 2003; 
Yan, Madgwick & White 2003; Kravtsov ct al. 2003). Note 
that this assumption is non-trivial, as it implies that galaxy 
formation is largely a local process in individual dark mat- 
ter haloes. This assumption was also made in many of the 
semi-analytic models based on Monte-Carlo halo merging 
trees (e.g. Kauffmann, White & Guidcrdoni 1993; Cole ct 
al. 1994; Sommerville et al. 1999), because in such models 
the merging histories are statistically the same for all haloes 
of the same mass at a given time Berlind ct al. (2003) found 
that this assumption is consistent with their numerical sim- 
ulations. 

To test the validity of this assumption, we use the halo 
occupation model recently developed by Yang et al. (2003a), 
which is based on the conditional luminosity function (here- 
after CLF) of galaxies in dark haloes of given mass, to 
predict how the galaxy luminosity function (hereafter LF) 
changes with large-scale environment. In this model, the 
change of the LF with large-scale environment is entirely 
due to the change of the halo mass function with large-scale 
environment. Since the conditional mass function of dark 
matter haloes can be accurately obtained from cosmological 
AT-body simulation, the CLF model can be used to make 
accurate predictions for how the galaxy LF changes with 
large-scale density field. We also show that, when averaged 
over a large volume, the predicted galaxy luminosity den- 
sity is closely correlated with the underlying mass density. 
Therefore, the luminosity density can be used to rank the 
mass density, and the model predictions presented here can 
be tested using large galaxy rcdsliift surveys, such as the 
two-degree Field Galaxy Redshift Survey (2dFGRS; Colless 
et al. 2001) and the Sloan Digital Sky Survey (SDSS; York 
et al. 2000). 

The paper is arranged as follows. In Section 2 we re- 
view the CLF model and describe how it can be used to 
calculate the LF of galaxies in different environments. Sec- 
tion 3 presents the model predictions for how the galaxy 
LF changes with large-scale density field. In Section 4 we 
use our model to study the correlation between the galaxy 
luminosity density and the underlying mass density, and dis- 
cuss how the predictions of our CLF formalism can be tested 
with observations. We summarize our results in Section 5. 



2 FROM THE CONDITIONAL LUMINOSITY 
FUNCTION TO GALAXY LUMINOSITY 
FUNCTION 

In Yang, Mo & van den Bosch (2003a, hereafter Paper I), 
we developed a formalism, based on the conditional lumi- 
nosity function $(L[M), to link the distribution of galaxies 
to that of dark matter haloes. We introduced a parame- 
terized form for ${L\M) which we constrained using the 
LF and the correlation lengths as function of luminosity. 
In van den Bosch, Yang & Mo (2003a, hereafter Paper II), 
we extended our model by constructing separate CLFs for 
the early- and late-type galEixies. The CLF formalism de- 
veloped in these two papers has subsequently been used to 
investigate galaxy clustering as function of luminosity and 



type (Yang et al. 2003b), to put constraints on cosmological 
parameters (van den Bosch, Mo & Yang 2003b), to con- 
strain redshift-evolution in the halo model (Yan, Madgwick 
& White 2003) , and to characterize the population of satel- 
lite galaxies (van den Bosch ct al. 2004). In this paper, we use 
the CLF formalism to predict how the luminosity function 
of galaxies depends on large-scale environment. For com- 
pleteness, we briefly summarize in the following the main 
ingredients of the CLF formalism, and we refer the reader 
to papers 1 and II for more details. 

The conditional luminosity function ^{L\M)dL gives 
the average number of galaxies with luminosities in the range 
L±dL/2 that reside in haloes of mass M. It is parameterized 
by a Schccliter function: 



$(L|M)dL = I- f -^Y exp(-L/L*)dL, 



(1) 



where L* = L*{M), a = ce{M) and $* = $*(M) are all 
functions of halo mass M. Following Papers I and II, we 
write the average total mass-to-light ratio of a halo with 
mass M as 



(2) 



which has four free parameters: a characteristic mass Mi, 
for which the mass-to- light ratio is equal to {M/L)o, and 
two slopes, 7i and 72, that specify the behavior of {M/L) 
at the low and high mass ends, respectively. A similar pa- 
rameterization is adopted for the characteristic luminosity 
L*{M): 



M 



L*{M) 



with 



r(a + 2) 
r(Q-f 1,1) 



(3) 



(4) 



Here r(x) is the Gamma function and r(a, x) the incomplete 
Gamma function. This parameterization has two additional 
free parameters: a characteristic mass M2 and a power-law 
slope 73. For a(M) we adopt a simple linear function of 
log(M), 



a(M) = ais+V log(Mi5), 



(5) 



with M15 the halo mass in units of 10^^ /i~^Mq, ais — 
q(Mi5 = 1), and 77 describes the change of the faint-end 
slope Q with halo mass. Note that once a{M) and L* (M) 
are given, the normalization of the conditional LF, 3'*(M), 
is obtained through equations (1) and (2), using the fact 
that the total (average) luminosity in a halo of mass M is 



(L)(M) 



f 

Jo 



#(L|M) LdL = !>* L* r{a + 2). 



(6) 



Finally, we introduce the mass scale Mmin below which we 
set the CLF to zero; i.e., we assume that no stars form inside 
haloes with M < Mmin- Motivated by reionization consider- 
ations (see Paper I for details) we adopt Mmin = 10^ Mq 
throughout. 

In order to split the galaxy population in early and 
late types, we follow Paper II and introduce the function 
/iate(i, M), which specifies the fraction of galaxies with lu- 
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minosity L in haloes of mass M that are late-type. The CLFs 
of late- and early-type galaxies are then given by 



^ute{L\M)dL : 
and 



/late(i,M)*(L|M)dL 



(7) 



Nearly (L|M) dL = [1 - /iate(i^, M)] *(L|M) dL . 



(8) 

As with the CLF for the entire population of galaxies, 
$iatc(i|M) and $eariy(I'|M) are constrained by 2dFGRS 
measurements of the LFs and the correlation lengths as 
function of luminosity. We assume that fig,te{L,M) has a 
quasi-separable form 



/iata(i^, M) = g{L) h{M) q{L, M). 
Here 

1 if g{L) h{M) < 1 



q{L,M) 



if g{L) h{M) > 1 



g{L)h{M) 

is to ensure that /iate(i, M) < 1. We adopt 
j'iato(L) /p°°<i>(L|M)n(M)dM 



9{L) = 



${L) ^{L\M) hiM) n{M) dM 



(9) 



(10) 



(11) 



where n{M) is the halo mass function (Sheth & Tormen 

1999; Sheth, Mo & Tormen 2001), l>iato(L) and <I(L) corre- 
spond to the observed LFs of the late-type and entire galaxy 
samples, respectively, and 



h{M) — max ( 0, min 



log(M/M„) 

log(Mi/Ma) 



(12) 



with Ma and two additional free parameters, defined 
as the masses at which h{M) takes on the values and 
1, respectively. As shown in Paper II, this parameteriza- 
tion allows the population of galaxies to be split in early- 
and late-types such that their respective LFs and clustering 
properties are well fitted. 

In Papers I and II we presented a number of difi^erent 
CLFs for different cosmologies and different assumptions re- 
garding the free parameters. In what follows we focus on 
the flat ACDM cosmology with Q.m ~ 0.3, flA ~ 0.7 and 
h = Ho/{WO kms"^ Mpc"^) = 0.7 and with initial density 
fluctuations described by a scale-invariant power spectrum 
with normalization ug — 0.9. These cosmological parameters 
arc in good agreement with a wide range of observations, 
including the recent WMAP results (Spergel et al. 2003), 
and in what follows we refer to it as the "concordance" 
cosmology. Finally, we adopt the CLF with the following 



parameters: Mi = lO^'^^/i" 
Ma = 10"-2'*ft-i Mo, Mb -■ 
124/1 (M/L)o, 71 = 2.02, 72 



M.:., M2 = 10'^ 
= 0.30, 73 = 0.72, 



"'^h-' Mo, 
(M/L)o = 
r? = -0.22 



and ai5 = —1.10. As shown in paper II, this model (referred 
to as model D) yields excellent fits to the observed LFs and 

the observed correlation lengths as function of both lumi- 
nosity and type. We emphasize, however, that our results 
are not sensitive to uncertainties in the CLF; had we chosen 
models A, B or C in paper II, instead of model D, the results 
presented below would have been virtually identical. Since 
$(L|M) is the average number of galaxies per unit luminos- 
ity in a halo with mass M, one can obtain the galaxy LF, 
$(L), from the halo mass function, n{M), according to 



$(L) 



-I 



^{L\M)n{M) dM. 



(13) 
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Figure 1. The contribution to the total luminosity function (solid 
curve) by haloes in various mass ranges: M/[h^-^ Mq] < 5 x 10^*^, 
5 X 10" - 10", 10" - IO12, 10" - 1013, IQi-^ - 1014, 10" - IQis, 
> IQI^ (broader curve corresponds to larger mass). 



Fig.l shows how haloes of different masses contribute to 
the total luminosity function. Note that the shape of the 
CLF changes significantly with halo mass, and that the faint 
end of the LF is dominated by galaxies hosted by low-mass 
haloes. 

If we make the assumption that the CLF is statistically 

independent of the large scale environment, the galaxy LF 
in a region of overdcnsity 5 follows from 



$(L|5) = / $(L|M)n(M|5)dM. 



(14) 



Here n{M\S) is the conditional mass function of dark mat- 
ter haloes, which gives the number density of haloes as a 
function of halo mass in an environment with average mass 
overdensity S = [Ai — A4]/M (with M the total mass in 
volume V, and the mean mass in all volumes V). Thus, 
in the CLF formalism, the (5-dependence of the galaxy LF 
enters only though the conditional mass function n{M\6). 

Fig. 2 shows the conditional mass functions for haloes 
in several representative environments. These functions are 
derived from a high-resolution A'^-body simulation which fol- 
lows the motions of 512^ particles with a P'^M code in a 
100/i~^Mpc box (see Jing & Suto 2002 for details), assuming 
the ACDM 'concordance' cosmology specified above. Dark 
matter haloes are identified using the standard friends-of- 
friends algorithm with a linking length of 0.2 times the mean 
inter-particle separation. As one can see, the shape of the 
conditional mass function is quite independent of 5 at the 
low mass end, while the shape at the high-mass end depends 
significantly on S. The break in the conditional mass func- 
tion at the high-mass end occurs at a lower mass for lower 
5. This reflects the fact that the formation of massive haloes 
is suppressed in low density regions due to bias (e.g. Mo 
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log[M/(h-iMj] log[M/(h-'Mj] 

Figure 2. The conditional mass functions of dark matter lialoes in various environments, specified by the mass overdensity 5 in a 
spherical volume with radius R = 8h~^ Mpc (left-hand panel) and R = 4:h~^ Mpc (right-hand panel). 



& White 1996; Gottloeber et al. 2003). For a given S, the 
break occurs at smaller mass if the volume used to define 6 
is smaller. 

Ideally, if we have an accurate model for the conditional 
mass function, we can combine it with the conditional lu- 
minosity function to construct an analytical model for the 
(5-dependence of the galaxy LF. Unfortunately, the model 
based on peak-background splitting (e.g. Cole & Kaiser 
1989; Mo & White 1996) is not sufficiently accurate for our 
purpose, even if ellipsoidal collapse is taken into account 
(Sheth, Mo & Tormcn 2001; Sheth & Tormcn 2002). The rea- 
son for this is that, when the total mass contained in a vol- 
ume becomes comparable to the mass of individual halos in 
consideration, as is the case for massive halos in low-density 
volumes, the peak-background splitting becomes inaccurate 
(e.g. Sheth & Tormcn 2002). We therefore decided to use 
the simulated samples constructed in Yang et al. (2003b) 
for our purpose. These simulated samples were obtained by 
populating dark matter haloes in A'^-body simulations with 
galaxies according to the CLF model described above. Here 
we use the sample constructed from one of the 100/i^^Mpc 
simulations, which is complete down to a &j-band luminosity 
of Lb. ~ W^h~^ L© (see Yang et al. 2003b for details). 



3 GALAXY LUMINOSITY FUNCTION IN 
DIFFERENT ENVIRONMENTS 

Fig. 3 shows the LFs of galaxies in regions of different 
mass overdensities S (defined over spheres with radius R = 
Sh^^ Mpc). As one can see, the shape of the LF at the faint 
end is quite independent of 6. This is due to the fact that 
most faint galaxies reside in low mass haloes, for which the 
shape of the conditional mass function n(A4\S) is only weakly 
dependent on S (cf. Fig. 2). The most pronounced difference 
between LFs in different mass density environments is the 
break at the bright end, which occurs at fainter luminosities 



CJ 



o 




9 10 
log[L/(h-2L3)] 



Figure 3. Galaxy luminosity functions in regions with different 
mass overdensities: —1.0 < <5 < —0.5 (open triangles), —0.5 < 
<5 < 0.5 (solid squares), 0.5 < <5 < 1.5 (asterisks), 1.5 < 5 < 4.0 
(solid circles), and <5 > 4.0 (open circles). The solid and dashed 
curves correspond to the best-fit Schcchtcr functions, fit over the 



luminosity ranges Li,j 
respectively. 



> lO^/i"^ Lq and Lf, > lO'-^/i"^ Lq, 



in lower-density regions. This, in turn, is simply a reflection 
of the fact that most bright galaxies reside in meissive haloes. 

Our CLF model also predicts a pronounced difference 
for the environment dependence of the LFs of early- and 
late-type galaxies. As shown in Fig. 4, the shape of the 
LF of early-type galaxies depends rather strongly on 5. For 
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the late-type galeixies, however, no such pronounced shape- 
dependence is predicted. Once again, this behaviour is easy 
to understand from the conditional halo mass frmction and 
the CLF: many of the fainter carly-typcs actually reside in 
clusters (see Paper II), whose abundance depends strongly 
on S. The majority of the faint late-type galaxies, on the 
other hand, reside in relatively low mass haloes, for which 
the shape of the conditional halo mass function is roughly 
independent of environment (cf. Fig.2). 

To further quantify how the LF depends on the mean 
density of the environment, we fit each of the luminos- 
ity functions by a Schcchtor form. For most of our discus- 
sion, we only fit over the luminosity range L > W^h^^ Lq 
{Mb- — Slog/i —17.2). In real redshift surveys, such as 
the 2dFGRS, galaxies fainter than this are observed only 
within a relatively small local volume, making it difficult to 
study their large-scale environment dependencies. As shown 
in Figs. 3 and 4, the Schechter function provides in general 
a reasonable fit to the LF over a large range in luminosity. 
This is not a trivial result, as the shape of the CLF varies 
significantly with halo mass and the conditional mass func- 
tion at the high-mass end changes significantly with 5. 

The left panel of Fig. 5 shows the characteristic lumi- 
nosity L* and the faint end slope a as functions of 5. In all 
cases, the characteristic luminosity increases with 5, reflect- 
ing the fact that massive haloes, which host bright galaxies, 
are preferentially located in high-density regions. The in- 
crease is only moderate: for the total sample it is about one 
magnitude from the lowest to the highest densities probed. 
The increase is more significant for early-type galaxies than 
for late-type galEixies. For late-type galax;ies, the increase 
becomes insignificant for 5 ^ 1. 

For the total population the faint-end slope a is roughly 
a constant, with a value about —1.1. The faint-end slope 
is also quite independent of 5 for late-type galaxies, with 
a ~ —1.25. In contrast, the faint-end slope for early-type 



galaxies changes from a ~ —0.3 for low-density regions to 
a ~ —0.8 for 5 ^ 10. As mentioned above, this is mainly 
due to the fact that a large number of relatively faint early 
types reside in massive haloes whose mass function depends 

strongly on 5. 

If we fit the luminosity function over the entire range 
with L ^ IQ'-^h-'^l.Q {Mb- - 5 log ft ^ -13.5), the faint-end 
slope becomes steeper for low-density regions (cf. open and 
solid circles in loft-panel of Fig. 5. Note that in this case, 
the Schechter function is no longer a good fit (see Fig. 3) 
for low-density regions; the faint-end slope is much steeper 
than that of the extrapolation from the luminosity function 
at the brighter end. This departure starts at L ~ lO®/i~^L0 
and is more significant for early-typo galaxies. 

We have also made similar analyses for spherical vol- 
umes with a radius of 4 /!.~^Mpc, the results of which are 
shown in the right-hand panel of Fig. 5. A noticeable dif- 
ference with using volumes of 8 /i~^Mpc radius is that the 
characteristic luminosity increases faster with 5. The reason 
is that the conditional mass function at the high-mass end 
reveals a stronger dependence on 5 for smaller radius R (see 
Fig. 2). 



4 PREDICTIONS FOR OBSERVATIONAL 
TESTS 

So far we have investigated how the galaxy LF depends 
on the mass density averaged over spherical volumes of 
8/i~^ Mpc radius. Unfortunately, it is not easy observation- 
ally to obtain accurate measures for the mass density field 
in the Universe. Therefore, in order to facilitate an obser- 
vational test of the predictions presented above, we need to 
find a quantity that is (i) easy to derive from observations, 
and (ii) that can be used to rank the mass density. 

When averaged within a large volume, the number den- 
sity, 5g of galaxies (with luminosities above some value) and 
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Figure 5. The characteristic luminosity and faint-end slope a as a function of S, defined as the mean overdensity within individual 
spheres of radius i? = 8 /i~^Mpc (left panel) and ii = 4 ?i~^Mpc (right panel). Symbols connected with lines are results of the fit 
over the luminosity range L > lQ^h~^ Lq for the total (solid line), early-type (dashed line) and late-type (dotted line) samples. For 
comparison, the open circles in the left-hand panel show the results for the total sample when fitting over the more extended luminosity 
range L > 10^-5/i~2 Lq. 



the luminosity density, 5l, are both expected to be corre- 
lated with the underlying mass density. Such relations can 
be derived directly from the simulated catalogs constructed 
with the CLF (Yang et al. 2003b). Fig. 6 shows the results. 
As one can see, there is a tight correlation between the lu- 
minosity overdensity Sl and mass overdensity S. Here 5l is 
calculated using all galaxies with L > Lmin = W^h^^ Lq, 
but it is not sensitive to this lower limit as long as Lmin ^ 
lO^/i~^L0. Thus the luminosity density can be used to rank 
the mass density. The use of galaxy number density is trick- 
ier, because the scatter is quite large for low-density vol- 
umes (for example, (5g, here defined by all galaxies with 
L > l(fh~^ Lq, becomes zero at 6 < —0.8 in the right 
panel of Fig. 6) and because it depends on the luminosity 
range used in defining the galaxy overdensity (5g (since the 
number density of galaxies is dominated by the faint ones). 

In order to show that Sl can be used to rank 5, we first 
fit the mean pL-p relation by the following function 

y = Ae-'-a^^e"*/", (15) 

where X — 1 + S,y = l + SL, and A, x* and /3 are fitting 
parameters. For high-density cells {x ^ x*), y oc x'^ , while 
for X <Si X*, y decreases exponentially with decreasing x. 
For = 8 h~^Mpc, we obtain A = 1.0203, x» = -0.3053 
and /3 — 0.8035. Wc then use this mean relation to convert 
the distribution function of the mass density into a distri- 
bution function for the luminosity density and compare it 
with the luminosity-density distribution function derived di- 
rectly from the simulation. The result is shown in Fig. 7. 
Clearly, the luminosity-density distribution function is re- 
covered quite accurately with the mean relation (15), sug- 
gesting that the luminosity density can be used to represent 
the mass density. 



1 I 1 1 1 1 1 1 1 1 1 1 1 1 1 r 

R=8 h-'Mpc 




log[l+5] 

Figure 7. The solid and dashed curves show the distribution 
functions of the mass overdensity 5 and of the luminosity den- 
sity Sl, respectively. The dotted curve shows the conversion of 
the distribution of mass density to the distribution of luminosity 
density using the mean Sl-5 relation shown in Fig. 6. 



Fig. 8 shows how L* and a change with the luminosity 
density. The results are similar to those shown in Fig. 5, 
except that the range of Sl is stretched at the low-density 
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-2-1012 -2-1012 
log[l+(5] log[l+(5] 

Figure 6. The left panel shows the relation between luminosity density {Si^, calculated using all galaxies resolved in the simulation) and 
mass density (S), while the right panel shows the galaxy number density {5g, calculated using galaxies with luminosity L > 10^h~'^ L©). 
All quantities are the averages within a sphere of radius R = 8 h~^Mpc, and normalized to the corresponding mean densities of the 
universe. The curves show the fits to the relations. Error bars for Sl and 5g represent the ranges of the quantities, while errorbars for 5 
are Icr standard deviations. Note that below a certain value of S, log(l + <5g) becomes — oo, while log(l + 5[,) is still finite. 
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Figure 8. The characteristic magnitude, M, , and the faint-end 
slope, a, of the best-fit Schcehtcr function as functions of ■ The 
results are shown for the total, early-type and late-type samples, 
as indicated. Schechter functions are fit to the luminosity func- 
tions over the range L > lO^ft."^ Lq. 

end relative to 6, because of the non-linear relationship be- 
tween Sl and S. As comparison, wc also show L* and a as 
functions of Sg, where 6g is based on the number of galaxies 
with L > W^h~^ liQ (Fig. 9). With such a luminosity hmit, 
the results based on Sg are very close to those based on Sl- 
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Figure 9. Same as Fig. 8, but as a function of Sg, the overdensity 
(in number) of galaxies with L > 10^h~^ Lq. 

5 DISCUSSION 

We have shown that the halo occupation model, which as- 
sumes that the luminosity distribution in a dark halo only 
depends on its mass, makes specific predictions about how 
the gala:xy LF changes with large-scale environment. The 
model predicts that the LFs for relatively bright galaxies can 
be fit reasonably well by a Schechter function, independent 
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of environment. In all cases the characteristic luminosity L* 
increases with local mass density. For late-type galaxies, the 
faint-end slope is quite independent of large-scale environ- 
ment while for early-types the value of a changes from ~ 
in low density regions to ~ —1 in high-density regions. The 
predicted LF for low-density regions shows significant depar- 
ture from the Schechter form at L ^ W^h~^LQ, and this de- 
parture is more significant for early-type galaxies. All these 
predictions have simple explanations based on the change 
of halo mass function with large-scale environment, and the 
change of halo occupation with halo mass. 

We have also shown that the luminosity density is 
tightly correlated with the mass density on large scales. 
Thus, our predictions can be checked using large galaxy red- 
shift surveys. At the moment, the study of the dependence 
of galaxy LF on local environment has mainly focused on 
comparing the LF of cluster galaxies with that of field galax- 
ies (see e.g. De Propris et al. 2003 and references therein). 
More related to the model predictions presented here are 
the analyses made by Hiitsi et al. (2002) and Bromley et 
al. (1998). Hiitsi et al. estimated galaxy LFs in spherical 
volumes with a radius 10 ft~^Mpc (in redshift space) within 
which the overdcnsitios of galaxy number are respectively 
5g < 0, < 5g < 1, and 5g > 1, and found that the faint- 
end slope a is about —1.1 for all the three cases, but that 
the characteristic luminosity brightens by about 0.4 magni- 
tude from the lowest-density sample to the highest-density 
sample. These results are consistent with our model predic- 
tion. Bromley et al. found that the faint-end slope of the LF 
of early-type galaxies (defined according to spectral type) 
depends strongly on local density, with a increasing from 
^ —0.4 in high density regions to ^ 0.2 in low density 
regions. This trend is also consistent with our model pre- 
diction. Unfortunately, current observational results are not 
yet sufficiently accurate to give a stringent constraint on 
the model. The situation will soon change. With the use of 
large redshift surveys of galaxies, such as the 2dFGRS and 
SDSS ^ , one cam estimate to high accuracy the galaxy lumi- 
nosity functions in regions of different luminosity densities 
or galaxy number densities (Croton et al. in preperation), 
facilitating an accurate comparison between model and ob- 
servations 

If our model prediction agrees in detail with observa- 
tional data, it will add strong support to the assumptions 
made in the theory: the 'concordance' cosmology (which de- 
termines the properties of the halo population) and that 
galaxy properties only depend on halo mass. If we assume 
that the halo population in the real universe is not very dif- 
ferent from that predicted by the 'concordance' cosmology, 
any significant difference between our model prediction and 
observation would have important implications for the for- 
mation of galaxies in the cosmic density field. For example, 
the prediction that the faint-end slope of the LF should not 
depend significantly on large-scale environment is due to the 
fact that most of the faint galaxies in our CLF model are in 



t During the final stages of tfiis project, Hoyle et al (2003) pub- 
lished a paper based on SDSS data that adresses the environment- 
dependence of the LF. Their results are qualitatively in good 
agreement with our predictions, though we delay a more detailed 
comparison to a future paper. 



low-mass haloes. A significant discrepancy between model 
prediction and observation would therefore require a CLF 
model in which many faint galaxies arc in massive haloes. 
This change of CLF is, however, strongly constrained by the 
observed weak clustering of faint galaxies. Thus, any sig- 
nificant discrepancy would strongly indicate the assumption 
that the galaxy properties only depend on halo mass is in- 
correct, invalidating the whole halo occupation approach. 
The other prediction of the model, that the characteristic 
luminosity increases mildly with local density, is due to the 
fact that in our CLF model bright galaxies are preferentially 
found in massive halos. Here any discrepancy with observa- 
tional data would therefore require a change in the fraction 
of bright galaxies in massive haloes. Such change will lead to 
a change in the correlation length as a function of luminos- 
ity, and is subject to sringent observational constraints. As 
mentioned above, the other CLF models we have obtained 
give very similar results, and so any significant discrepancy 
would mean that either galaxy properties depend signifi- 
cantly on other properties of halos in addition to the mass, or 
our parameterization of the CLF is not sufficiently general. 
As discussed in Paper I, our parameterization is motivated 
by physical considerations, and is quite general as long as 
the CLF is assumed to be a smooth function of halo mass. 
Thus, should a sigificant change in the CLF be required by 
the observation, it would have a significant impact on our 
understanding of galaxy formation in dark matter haloes. 

As far as an accurate comparison with observation is 
concerned, there arc several effects to be taken into account. 
The first is redshift distortion. Since bright galaxies in mas- 
sive clusters are expected to have large velocity dispersion, 
some of these galaxies may appear in 'low-density' regions in 
redshift space, weakening the dependence on 5. This effect 
was soon in our experiment where S was defined in redshift 
space. The second is concerned with the definition of over- 
density. In a magnitude-limited sample, volumes at large 
distances from us contain only bright galaxies, and so the es- 
timate of the overdensity in such a cell has to rely on a small 
number of bright galaxies. Since the shape of the luminosity 
function changes with S, the correction of the incompleteness 
due to the magnitude limit in individual volumes must take 
this effect into account. Finally, since the Schechter func- 
tion is not a perfect model, the fitting parameters depend 
on the luminosity range used for the fit. All these effects can 
be taken into account properly with the use of mock cata- 
logs by applying the same analyses to both the simulated 
and observed samples. We plan to present such a detailed 
comparison with observational data in a forthcoming paper. 

It is also interesting to compare the predictions pre- 
sented here, based on the CLF, with those based on semi- 
analytical models for galaxy formation (e.g. Kauffmann et 
al. 1999; Somerville & Primack 1999; Benson et al. 2002; 
Mathis et al. 2002). Mathis & White (2002) have shown 
that the characteristic luminosity of the galaxy luminosity 
function in their semi-analytical model increases progres- 
sively from low-density to high-density environments while 
the faint-end slope becomes slightly shallower. These results 
are in qualitative agreement with our predictions presented 
here. As already emphasized in paper II, this suggests that 
the halo occupation statistics as described by our CLF nicely 
fit within the standard framework of galaxy formation. 
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